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Based on a Perturbative QCD analysis of the deuteron form factor, a 
model for the reduced form factor is suggested. The numerical result is 
0^ | consistent with the data in the intermediate energy region. 
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S I. INTRDUCTION 

Exclusive processes involving the hadron at large momentum transfer were first studied 
in perturbative QCD(pQCD) by Brodsky and Lepage 0. Analysis of the deuteron form 
factor in the intermediate energy region |2| revealed that QCD could strongly affect the 
behavior of the deuteron electromagnetic form factors when Q 2 is the order of several GeV 2 . 
It was pointed out in Refs. || [| that the domain for leading-power pQCD predictions 
for the deuteron form factors is Q 2 3> 2MdAQCD ~ 0.8GeV 2 where a calculation with 
the Paris potential U shows explicit deviation, although it can explain the data well for 
Q 2 < lGeV 2 (See Fig. 1). In this domain the deuteron form factor can be written to the 
leading order in 1/Q 2 convolution: 

F d (Q 2 )= l\dx}[dy]^{y v Q)Tl q+ ^\x l) y j ,Q)^ d {x l) Q) , (1) 

J 



"Mailing address. 



1 



where the distribution amplitude $d(xi,Q) is defined as || 

rQ _ 

$d{xi,Q) = J [d 2 k ± }y 6q/d (xi,k ±i ) . (2) 
The notation [dx] and [<i 2 k_|_] is denoted by 

n n 

[dx] = 8(1 - J2xi) J[ dxi , (3) 



i=l i=l 

[d 2 kj S 16vr 3 5(E^)n^ • (4) 

i=l i=l 

However, the calculation of the normalized T# 9+7 to leading order in a s (Q 2 ) would 
require the evaluation of over 300,000 Feynman diagrams involving five gluons. Farrar et al. 
[|7| have done perturbative calculations on the helicity zero to zero deuteron form factor and 
found that it is much smaller than the deuteron form factor data at experimentally accessible 
momentum transfer. 

In order to make more detailed and experimentally accessible predictions, it was suggested 
in Ref. || to define a reduced nuclear form factor by removing the nucleon compositeness, 

UQ2) " ^wk ' (5) 

The argument for each of the nucleon form factors, F^, is Q 2 /4 since, in the limit of zero 
binding energy, each nucleon must change its momentum from P/2 to (P + q)/2. 

For the reduced form factor of the deuteron one finds the asymptotic scaling behavior || 

Q 2 UQ 2 ) - (inj^j , (6) 

which can be compared with the available data in the large Q 2 region, although this pre- 
diction is only for asymptotic momentum transfer. Equation reminds us the reduced 
form factor of deuteron may be derived in a way which is similar to the meson case in a 
perturbative QCD calculation. 

The aim of this paper is to build a model to calculate the reduced form factor of the 
deuteron in the intermediate energy region. The remainder of this paper is organized as 
follows: In Sec. II we analyze the reduced form factor, fd{Q 2 ), and the deuteron wave 
function. A QCD inspired model is built in Sec. III. In Sec. IV the numerical results for 
fd{Q 2 ) is given from our model. The final section is reserved for summary and discussion. 

II. THE REDUCED FORM FACTOR AND THE WAVE 
FUNCTION OF THE DEUTERON 

In the case of electron-deuteron elastic scattering, the standard Rosenbluth cross section 
is written (in the laboratory frame) as 

^ = (%)Mou[A(Q 2 ) + B(Q 2 )tan 2 ( 9 -)} , (7) 
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where A{Q 2 ) and B(Q 2 ) are determined by Gc, Gm and Gq 

A(Q 2 ) = Gl + ^Gl + l^Gl, 



B(Q 2 ) = -v(l + v)Gl (9) 



with 77 = Q 2 /AMj. The deuteron form factor F d (Q 2 ) is defined as F d (Q 2 ) = ^JA(Q 2 ). 

As mentioned in Sec. I, it is helpful to study the reduced form factor since the effects of 
nucleon compositeness (represented by F^) have been removed from it. Equation (|5|) means 
that the deuteron form factor F d (Q 2 ) can be factorized into two parts, and reduced form 
factor fd(Q 2 ) can be regarded as the form factor of a composite of two point-like nucleons. 
This factorization was obtained by assuming 

V d = ^ body x i> N x i> N (10) 

in a simple covariant model |[L1|| . ip^ is the nucleon wave function and 4u* Ay is the usual two- 
body wave function of the deuteron. The equation of motion for ^f d (xi, kj_j) in the light-cone 
frame(LCF) is given by 

6 1 2 2 

[M 2_]T ±i + mi ]* d (x i ,k ±i ) = [[dy][<Ph]V(x u k M]yj ,h3)My3>3^ • (11) 

i=l X i J 

The factorized form of the deuteron form factor can be got by substituting Eq. (|ll|) into 
Drell-Yan formula ||: 

6 

F d {Q 2 ) = £ e - I l dx ] [d 2 kuW d (x h k ±i + (S ia - Xi)q x )y d {x h k ±i ) , (12) 

a=l 

where qj_ is absorbed by the a-th quark, q = (0, q~, q±) and Q 2 = q^_. Noting that the gluon 
is a color octet in the SU(3) color group, the single-gluon exchange between two color-singlet 
nucleons is forbidden. Thus the real kernel calculation requires the inclusion of other com- 



ponents rather than two-nucleons. Brodsky and Ji (TTIj suggested a simple covariant model 
to incorporate the quark structure of the nucleon. The hard kernel at large Q 2 was assumed 
to be the perturbative amplitude for the six quarks to scatter from collinear to the initial 
two-nucleon configuration to collinear to the final two-nucleon configuration, where each 
nucleon has roughly equal momentum. They argued that the dominant configuration for 
this recombination is the quark-interchange plus one-gluon exchange between two nucleons. 
Thus, roughly speaking, we can divide the kernel into two parts. One represents the inter- 
change of quarks and the gluon exchange between two nucleons, which transfer about half 
of the transverse momentum of the virtual photon from the struck nucleon to the spectator 
nucleon. Another part is the inner evolution of two nucleons. The first part leads to the 
reduced form factor of the deuteron and the latter leads to the form factors of two nucleons 
together with the factorized wave function mentioned above. 
The body wave function in Eq. (|10D can be written as 



»?*(». U) = ^ri-^^] (13) 
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by using Brodsky- Huang-Lepage prescription [12| from a harmonic oscillator wave function 



A'exp[—^a 2 r 2 ], in the rest frame. As mentioned above, the reduced form factor can be 
obtained 

i2\ _ n / r^ir^.i^t. 



fd(Q ) = D J [dx][dy]<j> J d (x, Q)t H (x, y, Q)(j) d (y, Q) , (14) 
where D is a kinematic factor. The body distribution amplitude 4>d(x, Q) is defined by 

M*, Q)= j [dk ± ]^ ody (x,, k ±l ) (15) 

with xi + x 2 = 1 and k± 1 + kj_ 2 = 0. The kernel tn(x,y,Q) is dominated by the quark- 
interchange plus one-gluon-exchange diagrams. Since the binding energy of the deuteron is 
very small, the lowest Fock state(iViV configuration) is dominant. 

III. QCD INSPIRED MODEL 

For the deuteron case, the matrix elements of the electromagnetic current can be 
written in terms of three form factors as 

G» yx = (P'\'\J»\PX) 

= -{^(QV* • e[P" + P"*] + G 2 (Q 2 )[eV* • q - e'^e • q] 

-G 3 (Q 2 )e ■ qe'* • q(P^ + P'»)/(2M 2 )} (16) 

with Q 2 = —q 2 , q = P' — P, and e = e\, e' = ey are the initial and final polarization vectors, 
respectively. \PX) is an eigenstate of momentum P and helicity A. The Lorentz invariant 



form factors G{ are related to the charge, magnetic and quadrapole form factors [10 



G E = G x + -t]Gq , 
Gm = G 2 , 

G Q = G 1 -G 2 + (1 + V )G 3 . (17) 

Perturbative QCD predicts [|lj that the helicity-zero to zero matrix element Gq dominates 
helicity amplitude at large Q 2 for lepton scattering on the deuteron. In the standard LCF, 
defined by q + = 0, q y = and q x = Q, we have the following relations approximately: 



G 2 — 2C?i — - „_ 1 _ /ri — ; —Gqq , (It 



2 

2P+(2^+ 1) 

G 3 = (19) 



and 

Ge '■ Gm '■ Gq = (1 — ^77) : 2 : — 1 (20) 



while Q Aqcd- Also, Calson and Gross [[T3| have shown that the LCF helicity-flip ampli 



tudes G+ and G\_ are suppressed by factors of Aqcd/Q an d (^-qcd/Q) 2 , respectively. 
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Equation (|l4j) shows that the reduced form factor fd(Q 2 ) is determined by the body distri- 
bution amplitude <f) d (x, Q) and the kernel which can be obtained from the quark-interchange 
plus one-gluon-exchange diagrams between two nucleons. To represent these diagrams, a 
model can be built by introducing a vector boson(color singlet) with an effective mass M b 



14| . The picture of this model is described by Eq. (|14f) , which implies that the formulation 
is similar to that of meson form factor and t# can be computed in the one-boson exchange 
approximation, replacing the gluon by a massive vector boson and coupling constant g s by 
effective coupling constant g e $. We suggest that the effective mass M b can be determined by 
the empirical scaling law || 15]: 



O 2 

(1 + -^)f d {Q 2 ) = constant (21) 

with m 2 = M 2 = 0.28 GeV 2 . 

The hard scattering amplitude can be obtained by calculating the diagrams shown in 
Fig. 2 : 

4MV 2 1 

t H (x,y,Q)= TT^ 11 n : o > 22 

V ' xyQ 2 + M b 2 -{x-y) 2 M 2 xQ 2 + (± - (1 - x) 2 M 2 ) V ; 

where M is the deuteron mass, and we have taken the nucleon mass to be half of M. The 
kinematic factor D is 

D = Jl + ^n + ^n 2 . (23) 



IV. NUMERICAL RESULTS 

Substituting Eq. fl2"2"|) into Eq. (|14|), numerical analysis can be done with several param- 
eters, a, M b and a e ff = gl^/Aix, in the expression. From the empirical scaling law, M b 
is around 0.5 GeV. Instead of inputting the parameter a e s, we normalize the amplitude at 
Q 2 = 2.5 GeV 2 data point. The variation of the reduced form factor Q 2 fd(Q 2 ) vs Q 2 , with 
different values of M b , is displayed in Fig. 3. The effective coupling constant a e ff increases 
as M b become larger. By varying the parameter a in the wave function we can get different 
behaviors of Q 2 fd(Q 2 ) in the intermediate energy region, Q 2 > 1 GeV 2 . The corresponding 
effective coupling constant, which decreases as a become larger, is shown in Fig. 4. 

It is shown from fitting the data that M b = 0.5 GeV, a = 0.21 GeV and a e ff = 0.15. The 
result, comparing with the calculations with the Paris potential and the experimental data 
[|I~5| , [I~6|| , is shown in Fig. 1. Our results reveal that our model with the above parameters 
can explain the deuteron form factor well for Q 2 > 1 GeV 2 . 



V. SUMMARY AND DISCUSSION 

The scaling law of the reduced from factor suggest that the dominance of Gq begins 
at Q 2 ~ 1 GeV 2 . Thus one can calculate Gq to predict the reduced form factor in the 
intermediate energy region. 
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However, it is a very complicated problem to directly calculate Gq , since there are over 
300,000 diagrams and the evolution of the deuteron wave function leads to the dominance of 
hidden-color state contributions in the very large Q 2 region due to the gluon exchange in the 
kernel. In fact, Farrar, Huleihel and Zhang J7| found that hidden-color degrees of freedom in 
the deuteron wave function might be important in order to fit the data. In this paper we have 
tried to build a model to calculate the reduced form factor in the intermediate energy region, 
instead of doing a full QCD analysis. The point of this model is that the reduced form factor 
fd{Q 2 ) can be evaluated in a way similar to the meson form factor. It is determined by the 
body wave function 4>d{x, Q 2 ) and a kernel with a massive boson exchange. Our results show 
that our prediction can fit the data well for Q 2 > 1 GeV 2 . To fit the data we have chosen : 
(i) the effective gluon mass M& ~ 0.5 GeV, which is consistent with the empirical law, (ii) 
the parameter in the deuteron wave function a = 0.21 GeV, and (iii) the normalization at 
Q 2 = 2.5 GeV 2 data point, which corresponds to an effective coupling constant a e ff = 0.15. 

In addition, we restrict ourselves to calculate the reduced form factor only in the inter- 
mediate energy region. One can't expect that Gq dominates the helicity amplitude in the 
low Q 2 region, say, Q 2 < 1 GeV 2 . On the other hand this picture can't apply to the form 
factor at very large Q 2 since the hidden-color state contributions may be important in that 
region, where the full evolution of the six-quark wave function is involved. 

This model could be improved by taking into account the contributions of G+ and G~\__ 
in the low Q 2 region and the hidden-color contributions in the very large Q 2 region. We 
believe this model can be generalized to other light nuclei. 
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Figure Captions 



Fig. 1 Structure function A(Q 2 ) of the elastic erf-scattering from our model(the solid line) 
with M b = 0.5 GeV, a = 0.21 GeV, and the effective coupling constant a e ff = 0.15. The 
dashed line corresponds to the Paris potential calculation. 

Fig. 2 The hard scattering diagrams. 

Fig. 3 Comparison of the Q 2 fd(Q 2 ) data with our calculations by using the different effective 
mass of the vector boson, M b , while fixing a at 0.21 GeV and normalized at the Q 2 = 2.5 
GeV 2 data point. 

Fig. 4 Comparison of the Q 2 fd(Q 2 ) data with our calculations by using the different pa- 
rameter in the wave function, a, while fixing Mb at 0.5 GeV and normalized at the Q 2 = 2.5 
GeV 2 data point. 
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